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SIMULATIONS OF COSMIC CHEMICAL ENRICHMENT 
WITH HYPERNOVA 

Chiaki Kobayashi^ 



Abstract. We simulate cosmic chemical enrichment with a hydrody- 
namical model including supernova and hypernova feedback. We find 
that the majority of stars in present-day massive galaxies formed in 
much smaller galaxies at high redshifts, despite their late assembly 
times. The hypernova feedback drives galactic outflows efflciently in 
low mass galaxies, and these winds eject heavy elements into the inter- 
galactic medium. The ejected baryon fraction is larger for less massive 
galaxies, correlates well with stellar metallicity. The observed mass- 
metallicity relation is well reproduced as a result of the mass-dependent 
galactic winds. We also predict the cosmic supernova and gamma-ray 
burst rate histories. 

1 Introduction 



o 

While the evolution of the dark matter is reasonably well understood, the evolu- 
Q tion of the baryonic component is much less certain because of the complexity of 
, the relevant physical processes, such as star formation and feedback. With the 
■ commonly employed, schematic star formation criteria alone, the predicted star 
formation rates (SFRs) are higher than what is compatible with the observed lu- 
^ minosity density. Thus feedback mechanisms are in general invoked to reheat gas 
and suppress star formation. We include both supernova and hypernova feedback 
in our hydrodynamical model in this paper (see Kobayashi, Springel & White 
_C3 2006b, hereafter K06b, for the details). 

2 Hypernovae 

Supernovae eject not only thermal energy but also heavy elements into interstellar 
medium. Thus the star formation history is imprinted in the chemical abundances 
of stars and galaxies. Different types of supernovae, i.e.. Type II and la supernovae 
(SNe II and la), produce different heavy elements with different timescales. Re- 
cently, it is found that hypernovae (HNe), which have more than ten times larger 
explosion energy (£^51 > 10), produce a certain amount of iron. 
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Fig. 1. Cosmic star formation rates for no feedback case (dotted line) and for the SN 
(dashed line) and HN (solid line) feedback with the metal-dependent cooling. See K06b 
for the references of Rest-frame U-band, IR, and Ha observations. Rest-frame UV ob- 
servations are plotted with dust correction (Steidel et al. 1999). 



Kobayashi et al. (2006a) calculated the nucleosynthesis yields for wide ranges 
of metallicity and energy, based on the light curve and spectra fitting of individual 
supernovae. With these yields, the elemental abundance ratios are in good agree- 
ment with observations from oxygen to zinc. Especially, the observed abundance 
of Zn ([Zn/Fe] ^ 0) can be explained only by a large contribution of HNe. We 
adopt these yields and energies in our hydrodynamical model. 

3 Cosmological Simulations 

We simulate the evolution of gas and stellar systems and the chemical enrich- 
ment from SNe II, SNe la, and HNe, from the cosmological initial condition with 
Ho = 70 km s"^ Mpc'^, = 0.3, = 0.7, Qi, ^ 0.04, n = 1, and erg = 0.9. 
We use an SPH code GADGET-2 by Springel (2005), and introduce the metal- 
dependent cooling rates (Sutherland & Dopita 1993) and chemical enrichment 
scheme by Kobayashi (2004). The initial condition is calculated in a 10h~^ Mpc 
cubic box with periodic boundary conditions with A^dm — A"gas = 96'^. We adopt 
the Salpeter IMF, the star formation timescale c* = 0.1, and the feedback neigh- 
bors A^FB = 405. 

Figure [T] shows the cosmic SFRs that are directly measured from the ages of 
stellar particles for the low resolution simulations. Supernova feedback (dashed 
line) decreases the SFR from z ~ 3, but metal-dependent cooling increases as large 
as no feedback case. For comparison. Figure [T] also shows observational estimates 
of the cosmic SFR density at different epochs. Note that these determinations are 
derived from the observed luminosity densities, and this involves uncertainties from 
dust extinction and completeness, as well as from the IMF. If we include hypernova 
feedback (solid line), the SFR starts to be suppressed from z ^ 6 onwards, and 
is overall smaller by a factor of 3 at < z < 3. The resulting SFR is in broad 
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Fig. 2. The cosmic SFR for the different galaxy mass with the stellar mass of 10" (solid 
line), 10^° (dashed line), 10^ (dot-dashed line), and 10* M© (dot-dot-dot-dashed line). 
The dotted line shows the total. Galaxies are identified by FOF at z = (a) and at each 
redshift (b), respectively. 



agreement with the observations that show a peak of log(SFR/[M0yr~-'^]) ~ — 1 
at z ~ 4. Since the metallicity dependence of the cooling rate clearly appears for 
[Fe/H] > — 2, metal-dependent cooling does not play such a prominent role in the 
HN feedback case, and enriched gas can remain hot without forming stars. 

Without feedback, 25% of baryons turn into stars, which is too large compared 
with observational estimates. With SN feedback, the stellar fraction reduces to 
10 — 15%, which may be consistent with observation. Recently, the observational 
estimate has been reduced to less than 10% (Fukugita & Peebles 2004), which 
may require larger feedback. The larger energy ejection by HNe could provide a 
solution. The present mean gas metallicity is [Fe/H] ~ —1.3 for SN feedback, and 
is reduced to ~ —1.8 for HN feedback. The mean stellar metallicity is almost solar 
for the SN feedback, and becomes sub-solar for the HN feedback. 

When and where do stars form? To answer this question, we break up the 
cosmic SFR history according to galaxy mass. In Figure [2^, the galaxies have 
been identified by FOF at z = for this plot, and thus these SFRs correspond 
to the age distribution of stars in the galaxies. For all galaxy masses, the SFRs 
show a peak around z ~ 3 — 4, and the majority of stars are as old as ^ 10 Gyr. 
On the other hand, in Figure we identify galaxies at each redshift, which are 
comparable to the observations of high redshift galaxies. This shows that most 
stars have formed in low-mass galaxies with 1O^~^M0 at high redshift z > 3. 
10^°Mq galaxies exist at high redshift z < 5, but 10^"'^ Mq galaxies appear only 
after z ^ 2. From these two figures we conclude that most stars have formed 
in dwarf galaxies before they merge to massive galaxies in our simulation. As a 
result of the hierarchical clustering of dark matter halos, such old stars belong to 
massive galaxies at low redshifts. 

How are heavy elements ejected from galaxies to the IGM? In the simulation, 
we can trace the orbit of gas particles over time. Exploiting this, we define as wind 
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Fig. 3. The wind fraction against the total mass (a), and the stellar metallicity against 
the ejected metal fraction (b). 



particles those that are not in galaxies now, but have been in galaxies before. In 
this simulation, ^ 10% of baryons turn into stars, ~ 10% of the gas stays in 
galaxies (~ 8% is hot), and ~ 20% is ejected as galactic winds. The rest, half of 
the baryons, never accretes onto galaxies. 

When we follow the orbits of gas particles, we can also examine from which 
galaxies the wind gas particles are ejected. This allows a measurement of the 
ejected wind mass from each galaxy. In Figure[3K, a clear relation is found between 
the wind fraction and the total mass. Winds are efficiently ejected from small 
galaxies, with ~ 80% of accreted baryons being ejected from Mtot ~ lO^^M© 
galaxies. A similar relation is also found for the ejected metal fraction, i.e. the 
ratio between the wind metal mass to the total metal mass. It is interesting that 
the wind fraction and the ejected metal fraction correlate well with the stellar 
metallicity (Fig[3]D). Based on this finding, we conclude that the origin of the 
mass- metallicity relation can be explained with galactic winds. 

The metal enrichment timescale depends on the environment. Figure shows 
the evolution of oxygen abundances in the gas phase. In large galaxies, enrichment 
takes place so quickly that [0/H] reaches ~ — 1 at 2; 7, which is consistent with 
the sub-solar metallicities of the Lyman break galaxies (large errorbar at 2; = 3, 
Pettini et al. 2001). The low metallicities of DLA systems (errorbars, Prochaska et 
al. 2003) are also consistent with our galaxies, provided these systems are dwarf 
galaxies or the outskirts of massive galaxies. The low [C/H] of the IGM (box, 
Schaye et al. 2003) can be explained if the IGM is enriched only by SNe II and 
HNe. The average metallicity of the universe reaches [0/H] —2 and [Fe/H] 
~ —2.5 at 2; ~ 4 , but reaches the same values at z ~ 3 in the IGM. 

In galaxies, metallicity of the cold gas increases with galaxy mass (Figl3)D), 
which is comparable to observations with a large scatter (solid line at z = 0, 
Tremonti et al. 2004; z = 2, Erb et al. 2006). The central cold gas shows a 
relation between galaxy mass and metallicity with shallower slope than observed 
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Fig. 4. (a) Redshift evolution of oxygen abundances of gas. The points show individual 
galaxies with the size representing the size of galaxies. The solid, dashed, and dotted 
lines are for the mean, the largest galaxy, and the IGM, respectively. See the text for the 
observational data (errorbars). (b) Mean metallicities of cold gas (T < 10* K) within 
10 kpc, plotted against the total stellar mass at ^ = 0, 1,2, and 5. See the text for the 
observational data (solid lines). 



in emission-line galaxies. For the stellar population, the observed mass-metallicity 
relation is well reproduced (Kobayashi & Arimoto 1999; Pahre et al. 1998; Gallazzi 
et al. 2005), and originates in mass-dependent galactic winds. These relations are 
present since 2^5. 

4 Cosmic Supernova Rates 

We provide the cosmic supernova rates that is calculated with the cosmologi- 
cal chemodynamical simulations in Figure El Recently, large scale and deep 
surveys can give the SNe la rates as a function of redshift (jPahlen et al. 20041 



Barris & Tonry 2005). From the evolution of SN la rates, the delay time of SNe 
la have been constrained as a few Gyr (jMannucci et al. '2 006 ) . However, the the- 
ory of the SN la progenitor (i.e., a binary system of a white dwarf) is not simple 
and the evolution of the binary system should be taken into account. The de- 
lay time corresponds to the lifetime of the binary companion, and the shortest 
time of SNe la is much shorter as ~ 0.1 Gyr for the double-degenerate system 



dTutukov fc Yungelson 1994| ), - 0.3 Gyr for jMatteucci fc Recchi(2001)[ s model, 
and 0.5 Gyr for our single-degenerate system ( Kobayashi et al. 1998D . 

In addition, because the occurrence of SNe la depends on the metallicity in our 
SN la progenitor scenario, the SNe la rate depends on the environments, i.e., the 
star formation and chemical enrichment histories of host galaxies. Therefore, not 
with one zone model of the Universe but with the chemodynamical simulations, 
we should predict the cosmic supernova rates. In our simulations, the effects of 
hypernovae are included, and we can also provide the cosmic gamma-ray burst 
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Fig. 5. The cosmic rates for SNe II (long-dashed line), SNe la (solid line), and GRB 
(short-dashed line) per luminosity (SNu, upper panel) and volume (lower-panel). The 
dotted line shows the observation of GRB rates (Yonetoku et al. 2004) . The open triangle 
and open circle are for the observations of SNe II rates, and the other symbols are for SNe 
la. See Kobayashi & Nomoto (2006, in preparation) for the observational data sources. 



(GRB) rates assuming that the progenitors of GRB are massive (> 40Mq) and 
low- metal ([Fe/H] < — 1) hypernovae. Our predictions are in good agreement with 
observations. 
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